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We present a procedure that allows the construction of the metric perturbations and electromag- 
netic four-potential, for gravitational and electromagnetic perturbations produced by sources in 
Kerr spacetime. This may include, for example, the perturbations produced by a point particle 
or an extended object moving in orbit around a Kerr black hole. The construction is carried out 
in the frequency domain. Previously, Chrzanowski derived the vacuum metric perturbations and 
electromagnetic four-potential by applying a differential operator to a certain potential ^ . Here we 
construct $ for inhomogeneous perturbations, thereby allowing the application of Chrzanowski's 
method. We address this problem in two stages: First, for vacuum perturbations (i.e. pure gravita- 
tional or electromagnetic waves) , we construct the potential from the modes of the Weyl scalars i/jq 
or tpQ. Second, for perturbations produced by sources, we express 'I' in terms of the mode functions 
of the source, i.e. the energy-momentum tensor or the electromagnetic current vector Ja. 



I. INTRODUCTION 

The gravitational perturbations of Kerr Black holes 
(BHs) are fully described by the metric perturbation 
(MP) hap- The latter satisfies the linearized Einstein 
equation, which is a set of coupled, linear, partial dif- 
ferential equations. Teukolsky |^ showed, however, 
that the curvature Weyl scalars V'o and ■04 each satis- 
fies a decoupled field equation, the " master equation" . 
Furthermore, this decoupled equation may be separated, 
leading to ordinary differential equations for the radial 
and angular parts. This leads to a great simplification 
of the problem of determining the gravitational pertur- 
bations. 

The problem of electromagnetic perturbations over a 
Kerr background has a similar status. The Maxwell equa- 
tions form a set of coupled, linear, partial differential 
equations for the four-potential (or for the Maxwell 
field Fap). In this case, too, Teukolsky j^] j2| derived sep- 
arable, decoupled, equations for the two Maxwell scalars 
tpo and ip2- 

For several problems, e.g. the calculation of energy and 
angular-momentum outflux to infinity, knowledge of the 
Teukolsky variables (e.g. or ip2) is sufficient. However, 
there are problems for which one needs the full perturba- 
tion field (i.e. the MP hap in the gravitational case, and 
Aa - or alternatively the full tensor field Fap - in the 
electromagnetic case). This includes, for example, the 
calculation of gravitational or electromagnetic self force 
acting on a point-like particle orbiting a spinning BH. 

In principle, each of the Weyl scalars "00 and '04 con- 
tains the full information on the gravitational perturba- 
tion in vacuum ||^ (up to a few non-radiative degrees of 
freedoms, e.g. infinitesimal changes in the BH's mass and 
spin). Chrzanowski Q developed a procedure which al- 
lows the determination of the general homogeneous (i.e. 



vacuum) solution for the MP hap, by applying a certain 
differential operator to the homogeneous solutions for ip^ 
or ?/'4. It was later shown js), however, that this operator, 
when applied to a particular solution t/jq or ?/;4 yields a re- 
sult hafj which is a valid vacuum solution, but yet it rep- 
resents a physically different gravitational perturbation. 
Let us rephrase this in a more explicit manner: Con- 
sider a vacuum gravitational perturbation characterized 
by a particular function ^4. Then, there exists a certain 
function from which hap can be constructed by ap- 
plying Chrzanowski's differential operator. This function 
^I^ satisfies the same Teukolsky equation as the function 
p~'*'04 (where p is a certain quantity defined below), but 
yet ^ does not coincide with the quantity p~'^ip4 of the 
gravitational perturbation under consideration. 

The same situation occurs in the case of electromag- 
netic perturbations. Here, too, the full information about 
the (radiative part of the) electromagnetic perturbation 
is encoded in each of the Maxwell scalars ipo and ip2 
Chrzanowski's method [Q allows the determination of the 
general, homogeneous solution for A^ by applying a cer- 
tain differential operator to the homogeneous solutions 
for ipo or ip2. However, this procedure, when applied to 
a particular solution ipo or (^2, yields a vacuum solution 
Aa which represents a physically different electromag- 
netic perturbation [^. 

In view of the above, the problem of constructing the 
MP hap (the four-potential Aa) from ipQ or -04 {ipo or 
1P2) includes two stages: First, construct the potential 
vl^ from -00 or -04 {ipo or (P2), and second, construct hap 
(Aa) from ^E". The second part is well known — this is 
Chrzanowski's procedure The goal of the present pa- 
per is to address the first part, namely, the determination 
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of from ipo (or from ipo in the electromagnetic case).0 
This problem was recently addressed, for gravitational 
perturbations, by Lousto and Whiting (LW) ||] in the 
case of a Schwarzschild BH. Here we provide the solu- 
tion to this problem in the Kerr case (in the frequency 
domain) . 

We shall consider here two different physical situa- 
tions: (i) pure gravitational or electromagnetic waves 
(i.e. a vacuum perturbation in the entire spacetime), 
and (ii) gravitational (or electromagnetic) perturbations 
produced by a (charged) object orbiting the BH. The 
first problem is fairly simple, but the second one, that of 
perturbations with sources, is a bit more involved. The 
explicit solution for ^> in this case of inhomogeneous per- 
turbations, which is th e pri mary goal of this paper, is 
summarized in sections Vll^ (gravitational case) and IX 
(electromagnetic case) below. 

The equations relating the potential ^' to the rele- 
vant Teukolsky variables were derived by Wald g for a 
general algebraically special, vacuum, background space- 
time. The reduction of these equations to the Kerr case is 
given in Ref. Q for the electromagnetic case and in Ref. 

for the gravitational case. Our goal is to determine ^' 
by solving these equations. For either the gravitational 
or electromagnetic case, there are two such differential 
equations relating \1/ to the Teukolsky variables: a radial 
equation (i.e. one including r -derivatives), which relates 
^' to ipo or ifQ, and an angular equation (i.e. one includ- 
ing 6'-derivatives), relating ^ to -ipi or ip2. In Ref. |^ LW 
elaborated on the angular equation [Eq. (2.7) therein], 
and constructed its solution in the Schwarzschild case (for 
gravitational perturbations). Here we shall elaborate on 
the radial equation [Eq. (2.6) therein, or its electromag- 
netic counterpart] . This in fact turns out to be a simple 
ordinary differential equation, which is not difficult to 
solve even in the Kerr case. 

The MP hafj and four-potential Aa constructed via 
Chrzanowski's method are given in the ingoing (or out- 
going) radiation gauge Q] . Barack and Ori [gj recently in- 
vestigated the local asymptotic behavior of the radiation- 
gauge haf3 (either the ingoing or outgoing one) near a 
point particle, by locally integrating the equations defin- 
ing this gauge. They found that in this gauge /iq^ cannot 
be well defied all around the particle; Instead, there is a 
line of singularity that emerges from the particle to either 
the ingoing or outgoing radial direction, over which haf3 
diverges. (This line forms a 1 -f 1 dimensional singular- 
ity set in spacetime.) One can choose to have a regular 



function hap at r > r particle, where r is the radial coor- 
dinate, but this will inevitably lead to a line singularity 
at r < r particle] and vice versa. (Barack and Ori demon- 
strated this in the simplest case, i.e. a static particle 
located at r = Tparticie in flat spacetime, but the same 
situation should occur also for moving particles in Kerr.) 
Although the analysis in Ref. Q was restricted to the 
gravitational case, it is easily extended to the electro- 
magnetic case as well. It shows that the radiation-gauge 
Aa also has a line singularity, either at r > rparticie or at 

^ ^ "^particle- 

The solution constructed here provides an independent 
demonstration to this pathology of the radiation-gauge 
quantities /iq,^ and Aa near a point source. Throughout 
this paper we shall assume that the source is confined 
to a range r^-^ < r < r^ax (consider e.g. a point mass 
moving on an elliptical or circular orbit). In both the 
electromagnetic and gravitational cases (and for either 
the ingoing or outgoing radiation gauge) , one may choose 
to integrate the equations governing \1/ from r > rmax to- 
wards smaller r values. Then ^ is perfectly regular at 
r > ''max; but it turns out that at r < r„iax, ^' is irreg- 
ular on an outgoing null ray emerging form the particle 
inwardly (i.e. in the past direction). Alternatively one 
may integrate these equations from r < rmin towards 
larger r values, in which case ^' is perfectly regular at 
T < fmin but at r > Tmin it dcvclops an irregularity on 
an outgoing null ray emerging form the particle. This is 
perfectly consistent with the above mentioned irregular- 
ity of the radiation-gauge hap and Aa , found earlier (for 
hap) by Barack and Ori Throughout this paper we 
shall mostly refer to the solution for '5 which is regular 
at r > Tniax but has a line singularity at r < rmax, which 
we shall denote ^'^ for concreteness. The analogous so- 
lution (which is regular at r < rmin but has a line 
singularity at r > Tmin) may be constructed in a fully 
analogo us ma nner, as briefly summarized at the end of 
sections VIII and IX. (In section ^ we briefly discuss 
the possible implications of this line singularity to the 
self-force problem.) 

The two solutions '5+ and yield two different so- 
lutions for the MP hap or the four-potential Aa, both 
for the "same" (e.g. the ingoing) radiation-gauge condi- 
tion, which we denote /i^^ , A'^ and /i^^ , A^ , correspond- 
ingly. To avoid confusion we emphasize that these two 
solutions (for either hap or Aa, and, say, in the ingoing 
gauge) essentially represent the same physical perturba- 
tion, and they differ by gauge. That is, the ingoing (or 
outgoing) radiation-gauge condition does not completely 
fix the gauge. LO 



*We shall restrict attention in this paper to the construction 
of hap or Aa in the ingoing radiation gauge from tpo or <^o, 
respectively. The analogous problem of constructing hap or 
Aa in the outgoing radiation gauge from 1/14 or ip2 may be 
treated in a similar way. 



nt appears, though, that the "-f" and "— " perturbations 
also differ by some (non-gauge) non-radiative component, i.e. 
the so-called "Z = 0" and/or '7 = 1" modes. 
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Since there is full analogy between the gravitational 
and electromagnetic cases, the detailed calculations pre- 
sented in the next seven sections will refer to the grav- 
itational case only. The electromagnetic perturbations 
may be treated exactly in the same manner. Only in sec- 
tion IX we shall return to the electromagnetic case and 



irregular on S 



present the final procedure of constructing for electro- 
magnetic perturbations. 

In section || we give the basic field equations and es- 
tablish some notation. Section III presents the basic cal- 



culation of 4* in the case of pure gravitational waves, 
expressing it in terms of the modes of V'o- This result is 
then further simplified in section IV, by expressing both 
ipo and ^' in terms of basis solutions (of the relevant ho- 
mogeneous Teukolsky equation) which admit a simple 
asymptotic behavior at either the large-r limit or at the 
event horizon (EH). In section ^ we develope the gen- 
eral homogeneous solution for Eq. ( |ll|) below, the above 
mentioned "radial equation" relating 5" to "00 j which is 
a 4th-order differential equation. This general homoge- 
neous solution is required later for the construction of the 
relevant inhomogeneous solution (for perturbations with 
sources) . 



In section VI we address the physical situation which 
provides the main motivation for this paper, namely, 
gravitational perturbations produced by sources (e.g. a 
point particle, or an extended object, in orbit around a 
Kerr BH). We construct the solution for "if (more specif- 
ically, ^'"'") in this case, using the general homogeneous 
solution constructed earlier in section In the first stage 
the potential is expressed in terms of the inhomoge- 
neous mode functions of tpo- Then we further simplify 
the solution, expressing ^ directly in terms of the mode 
functions of the energy-momentum distribution in space- 
time. (In both stages we also use the homogeneous basis 
modes of ip4 in our expression for ^P). Whereas our de- 
tailed construction refers to '1'+, the construction of 
proceeds in a fully analogous manner, and we provide the 
final result for '5^ as well. Some details of the calcula- 
tions are giv en in Appendices A and B. 

In section VII we study the domain of validity of the 
solution 5*+ (and similarly of ^'^). For a point particle 
we find that is regular everywhere except at a (1 -I- 1) 
surface spanned by outgoing principal null geodesies em- 
anating from the particle's worldline in the small-r (i.e. 
past) direction. We denote this surface S+. For an ex- 
tended object, the solution 5*+ (which is typically regu- 
lar throughout) is valid everywhere except in a domain 
denoted (again) S+, which is now a four-dimensional 
set (the definition of which is provided therein). In a 
fully analogous manner, the other solution 5*^ is, for 
an extended object, valid everywhere except in a four- 
dimensional set S^; and at the point-like limit E~ de- 
generates to a (1 -f 1) surface spanned by outgoing princi- 
pal null geodesies emanating from the particle's worldline 
in the large-r (i.e. future) direction, with becoming 



Section VIII provides a summary of the construction 



of 5*, in the gravitational case, for the benefit of the 
reader who wishes to implement this method in prac- 



tical calculations. Then, in section [X we return to the 



electromagnetic problem and summarize the procedure 
of constructing ^' in this case, leaving many details of 
the derivation to Appendix C. Finally in section ^ we 
give some concluding remarks. 



II. PRELIMINARIES 

Consider the spacetime of a Kerr BH with mass M 
and specific angular-momentum a. We shall use the 
standard Boyer-Lindquist coordinates (t, r, 6, ip), and fol- 
lowing Teukolsky we denote A = — 2Mr + a?, 
S = + a^cos^fl, and p = — l/(r — iacosO). The 
Newman-Penrose Weyl scalars tpo and ipi (correspond- 
ing to s = +2 and s = —2, respectively) satisfy two 
decoupled master equations. Defining 

we may formally write the two master equations as 

W±2 [^±2] = 47rET±2 , (1) 

where W±2 is the second-order partial differential opera- 
tor 

Ws = -A-'^a.iA^+ia,] + [{r^ + a^/A - sin^ 9] du 
+ AaMrA'^d^t + [a^ I A - sin^^ e]d^^ 
-iiur^ede[suvede\ -2s[a{r - M)/ A + i cos 6* sin" ^ 6*] 9^ 
-2s [M{r^ -a^)/A-r- iacose]dt + (s^ cot^ ~ s) , (2) 

and T±2 is the corresponding energy-momentum source 
term given explicitly in Refs. Teukolsky further 
showed that these two Weyl scalars may be decomposed 
as 



V'±2 = V RiTir)S^l"^{e)e'^"''^- 



Xmuj 



(3) 



where i?^™'^ and S^'^ are, respectively, solutions of tlm 
radial and angular Teukolsky equations (given below) .11 
The source terms are expanded in a similar manner: 



i{inip—ujt) 



(4) 



The radial and angular ordinary differential equations 
take the form 



*In case the spectrum is continuous, the sum over lo should 
be replaced by an integral. 
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p^r [^irM] =r^rM (5) 

and 

e^^-" =0. (6) 

Here P^™'^ and O^™'^ are second-order linear differential 
operators, given by 

pX,nu ^ A-'dr[A'+^dr] + V ^(r) (7) 

and 

e^™"' = sin-i6iae[sin6i(9e] + Vg{9) , 

where the potentials are 

Vr{r) = [{r^ + a?ftj- ~ AaMuimr + a^rn^ 

+2iams{r - M) - 2iMujs{r'^ - a^)]A~^ 
+2iujsr ~ a^io^ — A 

and 

Veiff) — c?uP' cos^O — / sin^ 6 — 2aujscos9 

-2mscos9/ sin^ 9 - cot^ 9 + s + A . 

The parameter A is Teukolsky's separation constant, 
which we write as 

A = A-s - |s| , 

where A is the separation constant used by Chan- 
drasekhar Q (often denoted JK there). The parameter A 
runs over all eigenvalues of the angular Teukolsky equa- 
tion (^. Throughout this paper we prefer to use the 
separation constant A rather than A, because the angu- 
lar equations for s ~ 2 and s — ~2 have the sarne set of 
eigenvalues A ^ (which is not the case for A). I Also, 
this will allow an easier connection with Chandrasckhar's 
formalism. 

Our goal is to construct the MP. In a vacuum spacetime 
(i.e. T±2 = 0), the MP in the ingoing radiation gauge 
can be derived from a potential '^irg 0, by applying 
to the latter a certain second-order differential operator 
[0. We shall formally denote this differential operator 
by TliRG, namely, 

hiBG = T^iRG [f/flc] , (8) 

where hmc denotes the MP in the ingoing radiation 
gauge (for brevity we shall often omit the spacetime in- 
dices of the MP). Similarly, the MP in the outgoing ra- 
diation gauge can be obtained from another potential 



^In the special case auj = 0, the separation constant A be- 
comes l{l + l) — s'^ + \s\; hence, common A also means common 
I. 



^OUT, through another differential operator Hqut @, 
namely 

hoRG = ^ORG [^org] ■ 

In this paper we shall only consider the case of ingoing 
radiation gauge, but the potential '^org (and hence the 
MP in the outgoing radiation gauge) may be constructed 
in a fully analogous manner. For brevity we shall use 
here the notation 4" = ^/_rg (this variable is denoted tpc 
in Ref. |||]), 11 = Uirg, and haf3 = hiRc, hence 

= n . (9) 

The function "if has to be a solution of the vacuum 
Teukolsky equation for '0_2 namely, 

W-2 [*] = . (10) 

In addition, it must satisfy the following differential equa- 
tion I § 

^0 = [*] , (11) 

where throughout this paper an over-bar denotes com- 
plex conjugation. Here D is the differential operator 

where l'^ is the standard outgoing Kinnersley's tetrad vec- 
tor (see e.g. Ref. [0). We use here the abbreviated no- 
tation D'^ = DDDD , and the same for other operators 
used below. 

Our goal in this paper is to construct the function ^> 
that satisfies equations (^0|) and (pi]). This will allow 
the construction of hafj, the MP in the ingoing radiation 
gauge, through Eq. (^). We shall first consider the case of 
pure vacuum gravitational waves in the entire spacetime. 
In this case we assume that ^pQ is known (it encodes the 
information on the gravitational waves). Subsequently 
we shall consider the case of gravitational perturbations 
produced by a point particle (or any other matter source) 
moving in the Kerr spacetime. In this case we shall as- 
sume that T^2i the s — +2 energy- momentum source 
term, is known. 

III. PURE GRAVITATIONAL WAVES 

We now consider the case of pure vacuum gravitational 
waves, namely, T±2 = in the entire spacetime. The 
information about the gravitational waves is given by 
means of the Weyl scalar ^pQ. Since we are dealing here 
with linear perturbations, it will be sufficient to consider 
a particular mode Xmu of ipo] The entire perturbation is 
then obtained by a superposition. Thus, we now assume 
that Tpo takes the decomposed form 
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and the radial function satisfies the vacuum radial equa- 
tion: 



(12) equation (||) with — m, — w, and s = —2; and correspond- 
ingly we may write it as cS'^'2"'"'~", where c is a constant 



P+^"' [RlTir)] = 0. 



(13) 



We shall now use the Teukolsky-Starobinsky relations to 
show that the desired solution of Eqs. (|l0|,[ri|) is 



*=p-lA2(i?t)4[A2^0], 



(14) 



where p is a constant to be determined later, and D'^ is 
the differential operator 



-dt + dr - (a/A)9^ 



For a single Xmco mode we define the " reduced" operators 
D.,nu. ^dr+ iK/A , Dl^ ^dr- tK/A , 

where 

K = am — (r^ + a^)ijj , 
such that for any functions f{r) and g{0), 



D 



f{r)9{0)e 



g^e)e'(m^~^t)D^^ [fir)] , 



and the same relation holds between the operators 
D\D}-^^. (Note that Dmuj and Dl^^^ are the same as 



the operators "Pq" and "2?J ", respectively, in Chan 
drasekhar's notation Using the decomposition 



Eq. (|lj) becomes 



(15) 



(16) 



and Eq. (|l^) now reduces to 

R+T{r) - {D,nu.)^[R^-T] ■ (17) 

Let us first verify that 4" [the complex conjugate of Eq. 
(H)] satisfies Eq. (|l0|). The Teukolsky-Starobinsky rela- 



tions (see e.g. ||9|]) imply that the radial function R^^{r) 



constructed in Eq. ( [L6D is a solution of the s = — 2 radial 
vacuum equation, namely 



(r)]=0. 



(18) 



The complex conjugate of the radial function i?^™'^ is 
a radial vacuum solution with the same A and s = —2, 
but with negative sign for m and lo, which we denote 
j^x,-m,-uj rpj^g angular Teukolsky equation is real, and 
(holding A fixed) is invariant under the simultaneous 
change of signs of s, m, uj. Hence, S'^™'^(6') is a real func- 
tion which is also a solution to the angular Teukolsky 



(whose value is unimportant to us). We find that 



^ = cR 



A , — m , 



\0)e 



-iimip — Lut) 



(19) 



Thus, 'f is indeed a solution to the s — —2 vacuum 
Teukolsky equation (10) - a solution characterized by 
the set of indices A, — m, —uj. 

We still need to check that _R:^™" constructed in Eq. 
( |l6| ) satisfies Eq. (^7|) [this would in turn imply that the 
expression ( [l^ ) satisfies Eq. (pi])], and to determine the 
constant p. In fact all we need to show is that 

(A„.)^ [A'iDlJ^iA^R'^rir)]] 

is a constant multiple of i?^™'^. This follows immediately 
by applying the two parts of Theorem 1 in Chapter 9 of 
Ref. [^. Consequently, there exists a constant p such 
that 

RlTir)=p-\D„,^)^[A'{Dlj^[A'RlTir)]] ■ (20) 

From the analysis therein it becomes obvious that p is 
the real constant 

p = A2(A + 2)2 - 8oj^X[a'^{5\ + 6) - 120^] 



(21) 

arrUMj. Note that the coefficient p de- 
pends on the mode. Ell 

As was mentioned in section |, the potential ^' must 
also satisfy an angular differential equation, i.e. Eq. (2.7) 
in Ref. Q . The compliance of the above constructed vac- 
uum solution with this additional equation is guaranteed 
by virtue of the following considerations: (i) According 
to the analysis in Ref. [|^, there must exist a solution 
to the three simultaneous equations [i.e. Eqs. ([lO|) and 
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and the angular equation]; and (ii) The solution 
15,[l6]) is the unique solution to the simultaneous equa- 
tions ( p^ ) and (11). For, any nontrivial solution to the 
homogeneous part of Eq. ( pi] ) will violate Eq. ( p^ (one 
can easily verify this, based on the general homogeneous 
solution to Eq. ( pH ) , constructed in section ^ below) . 



IV. FURTHER SIMPLIFICATION OF THE 
VACUUM SOLUTION 

Equations (|l5|,^6|) provide the full solution for ^. It 
is possible, however, to construct a simpler and more ex- 
plicit expression for the radial function R^'^ . The latter 



**This may easily be deduced 
A^(L>i„„)*A2 to both sides of Eq. (|2C 
drasekhar's theorem 1, as well as his! 
9 of Ref. [§). 

^^It is assumed that p is finite and non- vanishing for all real 



applying the operator 
) , and then using Chan- 
?q. (43) (both in Chapter 
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which is the vacuum Teukolsky equa- the basis solutions R 



satisfies Eq. 

tion for the s = —2 radial function i?^™". Since this is 
a second-order ordinary differential equation, its general 
solution may be spanned by any pair of independent so- 



\maj(in,out) 
±2 



are defined to be the ex- 



lutions. Let R 



Xrnijj{a) 
±2 



and i?^™'^^''-' be such two pairs of 



independent solutions (one pair for s — +2 and one for 
s = —2). Let H denote the operator which maps a vac- 
uum s — +2 radial function i?^™" into the corresponding 
function ^^^'^ of Eq. @, namely, 



act solutions of the corresponding radial equations, which 
satisfy the asymptotic form ( ^ , P5|) at the leading order 
in 1 /r (the same remark applies to the event-horizon ba- 
sis functions defined below.) 

One can easily verify that the operators D„iuj , do 
not mix positive and negative exponents or r». There- 



fore, 
forms 

H[R 



and H\R 



Xmuj{out) 



must take the simple 



\mcj{in) 1 
+2 I 



(26) 



For each mode Xmuj there must exist a constant 2x2 
matrix Ca such that 



H[R 



+2 



(r)] 



CijR_ 



4j) 



(22) 



where i,j run over the two basis states a and b. The 
problem thus reduces to the determination of the four 
constants Cy . 

There are two preferred bases, however, for which this 
matrix becomes diagonal and especially easy to calculate. 
One such basis is the pair of solutions characterized by 
positive and negative exponents of at large r. Here 
is a function of r, defined by 



dr,/dr={r^ +a'^)/A 



(23) 



(and given explicitly below). Note that —>■ oo as 
r — > 00. The other basis is that of positive and nega- 
tive exponents of r» at the EH (where r* — > —oo). These 
two bases are also preferable for the physical interpreta- 
tion of the solution, and for its construction via a Green 
function (described in section VI below). The asymp- 
totic behavior of the vacuum radial Teukolsky functions 
are given in e.g. Ref. |^ for all values of s, both at the 
limit of large r and at the EH. 

In what follows we shall describe the application of Eq. 
(p2|), and the determination of the required coefficients, 
for these two special bases. 



and 



rrr T-tXmcjiout 

H [R_^_2 



(27) 



and the problem reduces to the determination of the two 
constants C'™) and C^™*). These constants may be de- 
termined from the large-r asymptotic form of Eqs. ( p^ ) 
or (p7|). Ignoring terms of higher order in 1/r, we have 



= dr — ■ 



and A = r^. 

In principle, both Eqs. (16) and ( |l7| ) may be used for 
the determination of each of the coefScients (7(0"*)^ 
Notice, however, that when Dmuj acts on a,nd 
^mt^ on A^i?^™"^*"-*, the leading-order term propor- 
tional to (jj cancels out. Therefore, in these cases the 
operator effectively decreases the powers of r (by 1 at 
least), as dr differentiates this power of r. This leads 
to a complication, because then we cannot ignore the 
higher-order terms (in 1/r) in the operators Dm^j, Dj^^ 
and in A, and also the higher-order terms in the basis 



solutions R 



Xmuj{in,out) 
±2 



On the other hand, no such can- 



celation of the leading order term occurs when D 

-,Amaj(m) 



acts 



on R 



±2 



and Dt, 



on K_j_2 . instead, we get 



[R 



Xmuj{in) 1 
±2 I 



-2iujR. 



Xmcj{in) 
±2 



and 



A. Large-r asymptotic behavior 



^rnuj[R±2 I = 2tUjR^2 



Considering the large-r asymptotic behavior of the vac- 
uum radial functions R^'^ (r) and R^'^ (r) , we may take 
the two basic solutions (for each s) to be those of positive 
and negative exponents of r* . These two solutions take 
the asymptotic form 



R 



Xmuj{in) 
+2 



'/r 



R 



Xmuj{out) 
+2 



and 



X7nuj{out) 



(24) 



(25) 



(see and recall the factor p ^ oc r^ in the above def- 
inition of ip~2)- To avoid confusion we emphasize that 



It will therefore be convenient to calculate C^"^^ from Eq. 
(0) and C(°"*) from Eq. (|l|), by substituting in these 
equations 



(28) 



(with " a" standing for either " in" or " out" , as appropri- 
ate). Equation (|l7| ) then becomes 



and Eq. (|T^) reads 



loo; p r R_^2 
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Since Eqs. (^,^^ imply 
we obtain 

C(") = l/(16w'*) , C(™*) = 16uj'^/p. (29) 



B. Event- Horizon asymptotic behavior 



In a completely analogous manner, we can use for the 
expansion of i?^™'^ and R^'^ basis solutions character- 
ized by either negative or positive exponents of at the 
EH (the latter corresponds to — > — oo). These basis 
solutions take the asymptotic form 



R 



Xmuj(down) 
+2 



2—ikr* ryXmuj{up) ^ ikr* 



and 



Rl 



Xmuj{down) ^ ^2^— ikr* j^Xmuj{up) ^ ^ikr* 



(30) 



(31) 



(see 1^). Here k = uj — ma/{2Mr^), where r+ is the r 
value at the event horizon, given by 

r± = M ± (M^ - a2)i/2 . 

In this case, again, one can verify that the operators 
Dmuj , Dj^^ do not mix positive and negative exponents 
of r*. Therefore, 



and 



jj^j^Xmuj{up)-^ _ q(up) j^Xmujiup) 



(33) 



and the problem reduces to the determination of the two 
constants C^'*"'"") and C("p). 

The leading-order form of A = (r — ) (r — r_ ) is 

A = q6r , 

where 5r = r — and 

q = r+-r^= 2{AP - c^fl"^ . 

Correspondingly, the leading-order forms of -Dmc^ and 
^Xau near the EH are 

Djni^ = ^^^!± (^dr, - ik) , Dj^^ = + ' 

where we have used + = 2Mr^. However, when 
applying the operators , Dj^^ to the above basis so- 
lutions, it is most convenient to view r and in Eqs. 
( pO| , ^l] ) as two independent variables. In this context we 
have 



and 



D,nuj = dr + ^i^!± _ ik) 



Dl^^dr+'^-^{dr,+lk) 



The above basis functions all take the general form 
F{r)e^^''^* . For such functions we have 



-ikr^ 



where w = AkMr+ and a prime denotes d/dr. Note that 

whpTi n artsnn „„H Ot nn A 2 nAmw(<io«)n) 

wneii lyrnu acts un -n,^2 aiici j-^jyi^ on z_i ^-j_2 ; 

the leading-order term proportional to k cancels out, and 

we are left with higher-order terms that take the lead, 

which is an inconvenient situation. For this reason we 

shaU calculate C('*°'"") from Eq. (0) and C("p) from Eq. 

(|l6|). Using again the substitution ( p8| ) (this time with 

"a " standing for either "down" or "up"), Eqs. ( p^ and 

( |l6| ) become, respectively, 

j^Xmu)(down) _ Q(down) f^jj s^i^j^Xmuj{down)T^ ^-g^^ 

and 

For the first equation we need to calculate the quantity 

(i?™^)4[A2e-''="] . (35) 

A straightforward calculation yields (at the leading or- 
der) 

(i?„^)4[A2e-^^-"] ^ gA-2g-«fcr* ^ 

where 

Q — {w ^ 2iq){w + iq)w{w — iq) . 
Wc thus find 

^(down)^^IQ^ (36) 

For the second equation we need to calculate the quantity 

(i?L)'[A^e''"]- 

This is just the complex conjugate of expression (^sf), and 
we find 



(Dl„^)4[A2e*'="] ^ QA-^e''^" 



which leads to 



(37) 
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V. THE GENERAL HOMOGENEOUS SOLUTION 
TO THE 4TH-ORDER EQUATION 

The equation ( |ll| ) that determines the potential ^ is 
an inhomogeneous fourth-order hnear differential equa- 
tion. In the last section we constructed the relevant in- 
homogeneous solution of this equation in the case of pure 
vacuum perturbations, for each mode of the source term 
-00 [by "relevant" we refer here to the solution that also 
solves the vacuum Teukolsky equation (p^j. Later we 
shall also need the general solution of this fourth-order 
equation, in order to construct the inhomogeneous solu- 
tions relevant to non-vacuum perturbations. To this end, 
we shall now construct the general homogeneous solution 
to Eq. (0). 

This equation is in fact a trivial ordinary differential 
equation. Let us denote by ^ the null geodesies whose 
tangent is the null tetrad vector l'^ (namely, ^ are the 
members of the outgoing principal null congruence). Let 
7 be an affine parameter along the geodesies ^, namely 

Then for any function f{t, r, 6, ip), 

D[f] ^ I'^d.f = ^{0 . 



The differential equation (|11|) thus reads 



d'j 



4 (0 = ^0 



(38) 



Its general homogeneous solution is a third-order poly- 
nomial in 7, whose four coefHcients may be taken to be 
arbitrary functions of ^: 



^ = ^ bi{^)Y (homogeneous) . 



(39) 



We wish, however, to re-write this homogeneous so- 
lution more explicitly as a function of the four space- 
time coordinates. To this end we need to explicitly 
parametrize the null geodesies ^. From the definition 
of l^^, along each null geodesic ^ we have 



r* (r) 



Specifically we take 



-dr , u{r) = 



A 



r* (r) 



■ ln(r — T-i-) 



and 



i{r) = (a/q) ln[(r - r+)/(r - r_)] 



dr . 

■ ln(r — r_) 

(41) 

(42) 



The null geodesies ^ are thus parametrized by the three 
parameters 9q = 9, Iq = t — r^:{r), and ipo = ip ~ u{r), 
and the above general solution takes the explicit form 

3 

* = ^ bi{Oo, to, (/7o)r' (homogeneous) , (43) 

i=0 

where bi are arbitrary functions of their arguments. 

Later we shall also need the form of this general ho- 
mogeneous solution in the frequency domain. In order to 
comply with the decomposed form (|l^), for a particular 
mode Amw the arbitrary functions hi must take the form 



where Bi are four arbitrary constants (for each mode). 
Correspondingly the (homogeneous-solution) radial func- 



tion i?^™'^ is given by 



'{r) 



iu-<^r.) ^ ^.^z (homogeneous) 



i=0 



(44) 



One can easily verify that this solution indeed satisfies 
the homogeneous part of Eq. ([l7| ) , namely 

{D,nu>f[R'-2'^] = (homogeneous) . 

To this end, it is sufficient to note that for any function 

fir). 



Dr. 



f{r)e 



dr 



-i{rnu—u)r^ ) 



(45) 



dr ^ d6 ^ dt 

d'j ' dj ' d'j 



/A 



We choose the origin of 7 such that 7 = r along the 
geodesic. Then t, 9, ip along the geodesic are given by 



9o , t = to + r^{r) , (p = ipQ + u{r) , 



(40) 



where 9o,to,ip>Q are arbitrary constants, and r^{r) and 
u(r) are given by the two integrals 



VI. GRAVITATIONAL PERTURBATIONS 
PRODUCED BY SOURCES 

Consider now gravitational waves produced by a point- 
like particle that moves freely in a Kerr spacetime. For 
concreteness let us assume that the orbit is confined to 
the range r^in < r < r,„ax (but this assumption may 
be relaxed, at least partially, as we discuss in section^). 
The orbit needs not be equatorial. Of special importance 
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is the case of a circular orbit, rmhi = '''max = ''o- Alterna- 
tively, we may assume that the gravitational waves are 
produced by a finite-size mater distribution. In this case, 
too, we shall assume that the matter is confined to the 
range Tmin <r< r,^^y,. 

In the formalism used here the single function 5* is re- 
quired to satisfy two differential equations - Eq. (^0|) and 
the inhomogeneous equation These equations are 

mutually-consistent in the case the source term for Eq. 
( pi] ) is a vacuum gravitational field iIjq, but otherwise we 
should expect to have an over-determination. Therefore, 
in a spacetime with a matter source (either a finite-size 
or a point- like source) , we cannot expect to have a solu- 
tion to both Eqs. (|l(j ) and (11) in the entire spacetime, 
or even in the entire vacuum part of spacetime. 

The nonexistence of a global solution 5" can be demon- 
strated from another point of view. For a point particle 
in an otherwise-vacuum spacetime, Barack and Ori H] 
showed there is no global radiation-gauge solution hap 
around the particle. As was discussed in section |, one 
can construct a solution K^p which is entirely regular at 
f > f particle, but this solutiou wiU necessarily have a 
singularity in the range r < r particle, along a line em- 
anating from the particle. Alternatively one may con- 
struct a solution h~p which is regular in the entire domain 
T < T particle i but this solution wiU have a line singularity 
in the range r > r particle- 

Although it is not possible to construct a solution 
valid in the entire vacuum region, it is possible (for a 
point source) to construct a solution which is valid ev- 
erywhere throughout the vacuum region, except in a set 
of zero measure. (In fact there are two such solutions, 
those denoted and in section || above.) We shall 
now proceed to construct such a solution. Specifically we 
shall describe the construction of the solution ^1'+ (but 
the other solution '5^ may be constructed in a fully anal- 
ogous manner). This construction is applicable in both 
cases of a point source and an extended source [though 
in the latter case the domain in which 5*+ violates the 
required equations is no longer of zero measure, as we 
discuss in the next section.] 

As in the previous sections, we shall consider here the 
fimction sourced by a particular mode Xmuj of '00 ■ 
This function takes the decomposed form (15), and we 
need to construct the radial function i?^™'^. In the vac- 
uum region r > Tmax, -R^™'^ is just the solution described 



in Section IV. This was shown to be a valid solution of 
both Eqs. (|17|) and ( p^ (and this is the only vahd so- 
lution). Note that in this external region ip^ is made of 
outgoing modes only. 



(r) =A(™*)i?+2(r) (r>w) 



(46) 



where hereafter we denote the relevant basis functions for 
brevity as 



R 



±2 



R 



X'muj(out) 
±2 



, R 



±2 



R 



X'inLo{down) 
±2 



and similarly we use C"^ = C^^ 

^Xmuj{out) 



c- 



(jidown)^ (Only 



Rj^2 ^^"i R±2 '^1^^ tie relevant here, because 

these are the two homogeneous basis functions involved 
in the construction of the retarded Green's functions for 
the Teukolsky variables '0±2; see below). Therefore, in 
f > ''max the radial function of (for a particular mode 
Xmuj) is simply given by 



i?^"-(r) = C+^(°"*)i?l2(r) (r > w) 



(47) 



Consider next the extension of this solution into the 
range r < Tmax- Here R^'^{r) is not a vacuum solution 
(it fails to be a vacuum solution everywhere in r-aan < r < 
''max), and we can no kinger require R^'^ to satisfy both 



Eqs. ([TtI ) and (p^. E3 We therefore choose to extend 
jl^nuj jj-^^Q J, ^ Tniax as a solution of Eq. (|l^), and, for 
the time being, forget about Eq. (p8|). Nevertheless, in 
the next section we shall show that this procedure yields 
a valid solution ^>'^, which solves both required (time- 
domain) equations (^0|) and (|ll]), even at r < rmax - 
except in the domain (which is of zero measure for a 
point source). 

To construct the solution of Eq. ( |l7|) we proceed as 
follows. The source term in this linear differential equa- 
tion is R^'^{r), the radial function of -00- This function 
is in turn sourced by the energy-momentum distribution 
in spacetime. It can thus be expressed by means of the 
energy- momentum source term T^™" (r) , via the Green's- 
fimction method: 



R 



+2 



T 



+2 



'{r')G{r, r')dr' . 



(48) 



The Green's function G(r, r') is constructed from the two 
vacuum solutions admitting the desired asymptotic be- 
havior, namely, outgoing waves at large r and down-going 
waves at the EH: 



G(r, /) = A+{r')RX2{r)0{r - r') 

+A-{r')Rl2{r)e{r' -r), 



(49) 



where 9 denotes the standard step function, namely 
9{x) = for a; < and 9{x) — 1 for a: > 0. The functions 
yl+(r') and A^{r') are given by 



^±-i?^2/(AT4^[i?-2,i?+2]) 



(50) 



(all quantities evaluated at r'), where W [R j_2 , RX-2\ de- 
notes the Wronskian of the two homogeneous solutions 



*£For, extending R^^" into r < rmax as a solution oLEq. 
( p8| ) would automatically yield the external solution ( [l7| ) in 
this range too. But then (Dmi»;)*[A^™"] would necessarily be 
the analytically-extended vacuum function R^", which does 
not conform with the actual, non-vacuum, function R^^'^ in 

^min ^ f ^ '"max- 
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i?^2- Note that G, viewed as a function of r, is continu- 



ous at r = r', namely 

A+ir')R+,ir') = A-{r')R-,ir'). 



(51) 



The coefficient A*^""*' in the external solution is thus 
given by 



T^^'^{r')A+{r')dr' . 



(52) 



We now wish to construct a "Green- like function" 
H{r,r') such that the function R^'^ will be given by 



R-Tir)= / TlT{r')H{r,r')dr' 



(53) 



in analogy with Eq. (^). This will be a solution to Eq. 
( |l7| ) if H{r, r') satisfies the equation 



{Dr^u.nH{T,T')] = G{r,r') 



(54) 



(in which the operator Dmoj differentiates with respect 
to r, not r'). Motivated by the above form of G{r,r'), 
we assume a function H{r, r') of a similar form, 

H{r,r') = H+{r,r')e{r - r') + H-{r,r')e{r' ~ t) , (55) 

where H^(r,r') and H^(r,r') are smooth functions of 
their arguments. Equation ( p4[ ) is then satisfied if the 
following two conditions hold: (i) the two functions 
H^{r, r') satisfy 



(56) 



and (ii) H{r,r') is continuous and differentiable four 
times (with respect to r) at r = r'; In other words, the 
function 

y{r,r') = H+{r,r') - H^{r,r') , 

and its derivatives with respect to r up to 4th order, 
vanish at r = r': 



5" 2/ 



(r = /)=0, n = 0, . 



(57) 



(in which d'^y/dr^ = j/ is to be understood). Condition 
(i) guarantees the validity of E q. (|54| ) at r > r' and r < r' 
[the "-I-" and "— " cases in Eq. (|56[), respectively]. Condi- 
tion (ii) is required for the validity of Eq. ( [S^ ) at r = r'. 
[To see this, rewrite Eq. (|55|) as 



(r, r') = y(r, r')0{r - r') + H' (r, r') , 

and recall the continuity of G at r = r'.] 

The form of the general solution (r, r') to Eq. ( |5^ ) is 
obvious from the analysis in sections |^ and ^ above. We 



have a specific inhomogeneous solution C^A^{r')R%{r), 
and the general homogeneous solution (|4j) (in which the 
coefficients are now allowed to be arbitrary functions 
of r'); Hence the most general solution is 

3 

H^{r,r') ^ C^A^{r')R%{r)+e-'^"'''-'^''''^^Bf{r'y. 



1=0 



(58) 



It should be noted that since this is the most general so- 
lution to Eq. (17), this form must be satisfied by both 
and 5"^. The difference between these two solutions 
should emerge from the choice of the free functions Bf, 
which are to be determined by the boundary conditions. 
As we are considering here the solution 5'"^, the radial 
fimction _R^™"(r) must satisfy Eq. (E^) at r > Tmax- 
This is achieved by simply choosing Bf{r') = for all i, 
namely 

H+{r, r') = C+A+{r')Rt2{r) ■ (59) 
The internal part {r, r') has a more complicated form, 

3 

H-{r,r') = C-A-(r')i?l2(?-) +e^''(™-"''*)^Sr(r')r\ 

i=0 

(60) 

in which the four arbitrary functions B^{r') are to be 
determined by condition (ii) above, i.e. by matching to 
H'^{r, r') at r = r' . The function y[r, r') is given by 

y{ry) = C+A+{r')RU{r) - A' {r')R-_^{r) 

3 
i=0 

and we must impose Eq. (|5^). This might look problem- 
atic at first glance, because apparently the latter equation 
imposes five requirements on the four arbitrary functions 
B~[r'). However, one of these requirements is automati- 
cally satisfied. To see this, it will be convenient to rewrite 
Eq. d^) as 



0. 



0, 



(r 



(where (-Dmw )"[?/] = y is to be understood). Considering 
the case n = 4, the operator {DmujY annihilates the last 
term in the above expression for y{r, r'), and we have 

(A„J*[y] = A+ir')R+,{r) - A'{r')R-,{r) , 

which vanishes by virtue of Eq. (|5^). We can therefore 
re-express condition (ii) as 



{D„,^ny]=Q, 



n = 0, 3 



(r = /) . (61) 



The four arbitrary functions B^ (?*') should thus be deter- 
mined from the four conditions involved in this equation. 
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In Appendix A we show that these four functions may be 
expressed as 



Br{r') = C+A+{r') 
-C-A-{r') 



f,[r')R+^{r')+g,{r')—R\{r') 



where fi(r') and gi{r') are certain functions of r' exphc- 
itly specified therein. 



A. Further simplification of the solution 



and 

-a-{r') 



f,{r')RUr')+g.{r')—RU{r') 



f,(r')RZ,{r')+g.{r')^RZ,{r') 



Note that when expressed in this form the function 
H{r,r') — and hence also R^'^ir) — is invariant to a 
rescaling of i?^2 O'" ^-2 by constants. Therefore there is 
no need to require here a specific normahzation for these 
basis solutions. 



The procedure described so far for the construction 
of R/^'^(r) requires the two s = —2 homogeneous ba- 
sis functions -R^2i ^1^° the two s — +2 homoge- 
neous functions R^2- The latter functions are required 
for the determination of A^{r') (and their derivatives are 
involved in the Wronskian VF[i?^2: ^+2])- However, in 
principle ipQ can be determined from ^4 (and vice versa), 
and this implies that R^2 may be determined from R^2 ■ 
In Appendix B we undertake this goal and re-express 
A^{r') in terms of the functions R^2 ^-nd their first-order 
derivatives. We find 



A±(r') = (pC^W[RZ2.R- 



A{r')—RUr') + B{r')R^2ir') 



(62) 



where p is the parameter defi ned in Eq. (^l|), A and B 
are functions specified in Eq. (B7), and 



W[RZ2,Rt2] = const ■ A(r') 



(63) 



is the Wronskian of the two basis solutions R^2 (evalu- 
ated at r'). 

In the above construction of H{r, r'), A^ and only 
appear through their products A~^C~^ and A~C~. We 
therefore define 



a±(r') = C^A^ir'), 



and obtain 



a^{r')^ {pW[RZ2,Rt2]) ' 



A{r')—RT2ir') + Bir')RUr') 



(64) 



The functions H^{r,r') and Bj^{r') can now be re- 
expressed as 



H+{ry)=a+{T')Rl2{r). 



(65) 



F-(r, t') = a-{T')RZ2{r) + e-^("™*) ^ Biir'y , 



i=0 



(66) 



VII. DOMAIN OF VALIDITY OF THE 
CONSTRUCTED SOLUTION 

In the case of perturbations produced by matter 
sources, we must carefully examine in what parts of 
spacetime the construction of ^ and h^p is valid. First, 
the Chrzanowski's construction requires that the poten- 
tial 5" satisfies both equations (|lO|) and (pj]). These equa- 
tions are mutually consistent in vacuum, but are gener- 
ally inconsistent in the presence of matter. The matter 
source therefore leads to violation of at least one of these 
equations. This violation occurs not only in the region 
occupied by the matter, but also in certain vacuum parts 
of the spacetime .III 

At this stage it will be conceptually simpler to assume 
that the massive object that creates the perturbation has 
a finite size and a regular energy-momentum distribution 
7)11/(2;"). (The case of a point mass will then follow in 
a trivial manner.) Let us define E+ to be the collection 
of all points P in spacetime which have the following 
property: The null geodesic ^ (a member of the outgoing 
principal null congruence) passing through P intersects 
matter (more precisely, nonvanishing source term T+2) 
on its approach from P towards future null infinity. The 
collection of all other points of (the part r > r+ of) space- 
time is denoted In an analogous manner, we define 
E~ to be the collection of all points P for which in- 
tersects matter on its approach from P towards the EH, 
and S~ is the rest of (the part r > r+ of) spacetime. In 
a more pictorial language, imagine that light rays prop- 
agate along all geodesies ^ of the outgoing principal null 



^^To verify this, note that in the above construction of '1'+, 
the individual-mode radial function i?:^'J"(r) violates the ra- 
dial Teukolsky equation in the entire domain r < rmax — 
which in particular includes the vacuum domain r < Vmin- 
This violation follows from the non-vanishing of the coeffi- 
cients B~ . This does not necessarily mean that the time- 
domain Teukolsky equation ( |lO[ ) is violated everywhere in 
r < r-niax (in fact it does not, as we show below); but it does 
indicate the existence of a domain of violation that extends 
at any r value in r < Tmax- 
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congruence. Then S~ is the portion of spacetime "shad- 
owed" by the matter-energy distribution T+2(a;"), and 
is the non-shadowed part. Similarly, is the "past- 
shadowed" part, i.e. the part of spacetime that would 
be shadowed if the light rays were propagating along the 
null congruence from future to past (and from large r 
towards the EH), and E"*" is the rest of r > r+. By 
definition, S+ and E~ are pure vacuum domains. Note 
that E"*" contains the entirely- vacuum domain r > rmax, 
and generally it also extends through r < rmax into the 
other entirely-vacuum domain r < rmin — though the 
latter domain is not entirely contained in E"*". (Similarly, 
E~ contains the entire domain r < rmin, and generally 
extends through r > r^-m into r > Tmax)- Also, at the 
point-like limit E+ and E~ each degenerates to a (1 + 1)- 
dimcnsional surface, that emerges out of the particle's 
worldline in either the past or future (i.e. inward or out- 
ward) direction of ^. Hence in the point-like limit E+ or 
S~ cover the entire spacetime except a set of measure 
zero. On the other hand, when the object is extended, 
E+ and E" are four-dimensional sets. 

We shall now argue that the above constructed po- 
tential \E'+ — and the MP h'^„ constructed from it via 

Eq. (||) — are valid in the entire domain E+. We shall 
first establish this for the potential ^1'+, namely, we shall 
show that Eqs. (Iin) and (O) are satisfied throughout E+. 



Then we shall show that the MP solution h'^„ constructed 

ap 

from is valid throughout E+ too. (The same argu- 
ments apply to the validity of and h"^ throughout 

S-.) 



A. Validity of the constructed potential 

In the range r > rmax, the above construction of the 
"Green-like function" H{r,r') ensures that Eq. ( ^) i s 
satisfied by R^'^{r), and hence also Eqs. (0) and (ll8|). 
In the time domain this implies that satisfies Eqs. 
( p^ and (|l^) throughout r > rmax- 

In the range r < rmax we have constructed R'^'^{r) 
such that it satisfies Eq. ( p^ for all modes, hence in 
the time domain compliance with Eq. ( pT| ) is guaranteed. 
But Eq. ( p^ ) is violated throughout r < rmax [note that 
the homogeneous solution (H) violates Eq. (|l|)]. We 
shall now employ analyticity considerations to demon- 
strate that does satisfy the time-domain Teukolsky 
equation ( p^ throughout E+. 

Each mode Xmuj of is analytic throughout r > rmax 
by construction [cf. Eq. (47)]. Assuming convergence of 
the mode sum at r > rmax, wepjuay assume that 



itself is analytic at r > 



too 



Now, in the above 



construction we extended into the entire domain r < 
r-max as a solution of Eq. (pT|). This implies that 5*+ is 
analytic throughout E+, as we now show. 

Equation ( pT| ) is an ordinary differential equation along 
the null geodesies ^, which we may write as 



(67) 



We write its general solution explicitly (in a recursive 
manner) as ^+(7;^) = $1(7; with 



^niT,0 = c„(a + J <f>„+i(7';C)d7' n = 1...4 , 

ro 

(68) 

where $5 = -00 • Here Ci{^) [i — 1...4) are four arbitrary 
functions of the three variables 6'o,to,Vo- [Recall that 
the geodesies ^ are parametrized by the three quantities 
6'o,io,<i5o, defined through Eq. (^, that take constant 
values along the geodesic. Also recall that we have set 
7 — r] We take the lower integration limit to be, say, 
ro = 2rmax- The transformation from the coordinates 
{t, r, 6, (p) to (7, ^0, ^0, ^0) can be read off Eq. (^), and it 
is manifestly analytic everywhere in r > r+. Therefore, 
the analyticity of in the domain r > rmax implies it 
is analytic in (7, ^Oi^Oi Vo) as well. This in turn implies 
that all four functions Ci(^) = Ci{6o, to, ipo) are analytic in 
(^Oi toi^o)- Now, the function ^0 JS-jpresumably analytic 
everywhere in the vacuum region .til 



When the solution 
( |68| ) is restricted to the domain E+ , we observe that only 
vacuum points (7';^) are encountered in the integration, 
hence ^'0(7'; analytic. This immediately implies that 
given in Eq. (|6^ ) is analytic in (7; 6*0, to, fo) through- 
out E+, and hence also in (i, r, 9, ip). 

From the analyticity of (which implies the analyt- 
icity of *+) it follows that W-2 [^^] is analytic too. The 
vanishing of the latter at r > rmax therefore implies its 
vanishing throughout E+. We have thus established the 
compliance of with Eqs. (|l0| ) and ( |ll| ) throughout 
E+. 

It is easy to see why this argument fails at E+: The 
function tpo fails to be analytic at the point particle, or 
— in the case of an extended object — at the bound- 
ary of the matter distribution. As a consequence, along 



***For our argument to hold it is sufficient that the mode 



sum converges throughout some range r > r^ax > rmax , or 
even throughout some open interval of r values located some- 
where at r > rmax. 

^t^In the point-like case tpo is irregular at the particle's loca- 
tion. In the case of a smooth extended source, ipo will fail to 
be analytic at the boundary of the region occupied by mat- 
ter. But in both cases we may assume that ^po is analytic 
throughout the vacuum region. 
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each null geodesic ^ intersecting the source, will be 
analytic only up to the intersection point. Then will 
usually be non- analytic at the boundary of E+. There- 
fore we cannot expect to satisfy Eq. ( p^ in S+. The 
violation of the corresponding frequency-domain equa- 
tion (18) throughout r < r„iax indicates that Eq. (O) is 



indeed violated somewhere in T,'^ (and for any value of ', 
in this range). 

Finally we note that the compliance of 4*+ with the 
"angular equation" (i.e. Eq. (2.7) in Ref. Q) through- 
out E"*" may be deduced by exactly the same analyticity 
argument. 



and our goal is to construct the radial functions JZ^™". 
We shall now summarize this construction in the two 
different cases: (i) pure gravitational waves, and (ii) per- 
turbations with sources. 



A. Pure gravitational waves 

In this cases we assume that ipQ is given. This field is 
decomposed into modes too, 



B. Validity of the constructed metric perturbation 

Our goal here is to establish the validity of 
(constructed from via the Chrzanowski's method) 
throughout S"*" , despite the presence of matter in space- 
time. By "validity" we mean that (i) /i^^ satisfies the lin- 
earized vacuum Einstein equations, and (ii) the s — +2 
Weyl scalar constructed from it coincides with the orig- 
inal field V'o- To this end we use analyticity considera- 
tions, similar to those used above for analyzing the va- 
lidity of ^f^. Here we shall briefiy sketch these consider- 
ations, lol 



Consider first the validity of in 



r > r-max- To this 
end, expand (and -00) into modes. For a particular 
mode Xmuj, extend the r > rmax vacuum solution analyt- 
ically into r > Tmax- This extended solution represents 
a pure vacuum perturbation. Chrzanowski's construc- 
tion may now be applied to it, yielding the MP solution 
^+Ama; ^ Yl [^+Amt^j f^^. ^j^g mode under consideration. 



For each mode Xmu, R^'^{r) is a solution of the vacuum 
radial Teukolsky equation, and we assume this function 
is provided as a linear combination of two basis solutions. 
Two sets of convenient basis solutions are: (i) the large-r 
set , in which the basis solutions for i?^™'^ and i?^™'^ are 
given in Eqs. ( ^ , p5|) , and (ii) the EH set, in which the 
basis solutions for i?^™'^ ^^'^ given in Eqs. (^ , ^ ). 

Assume now that the information about -00 is given in 
terms of any two of the above four s = +2 basis functions, 
namely 

and the coefficients and A'^^^ are provided for each 
mode XmLo. (Here "a" and "6" denote either the large- 
r basis solutions, or the horizon basis solutions, or any 
combination of these two sets, e.g. ^^a"="out" and 
" 6" =" down" ) . Then, the corresponding radial functions 
il^nuj q£ ^ g^j.g simply given by 



Upon summation over the modes, we obtain a valid MP R^'^{r) = C'^'^A^'') R^^'^^°-\r) + C^''^ A^''^ R^_ 



4b) 



solution = H[\['+] in the range r > r^nax- 

Next consider the validity of the solution h'^^ = H [^f+J 

in the part r < Tmax of E+. From the analyticity of 5'+ 
throughout E+ (established above) it follows that ft,^^ is 
also analytic in this range. Recall also the analyticity 
of "00 throughout E+. The above criteria (i,ii) for the 
validity of a MP solution hap are both formulated in 
terms of analytic differential operators acting on hap. 
From the validity of these criteria in the range r > rmax 
it now follows that they must hold throughout E^. 



VIII. SUMMARY OF MAIN RESULTS: 
GRAVITATIONAL PERTURBATIONS 

Here we briefly summarize our procedure for construct- 
ing the potential for gravitational perturbations. We 
use the decomposition 



(r). 



The four coeflncients C^"'\ c'""*), C^^own) g^^^j (j{up) ^^^^ 
specifled in Eqs. (p9l|6l|3^). 



B. Gravitational perturbations produced by sources 

Here we consider the case in which the perturbation is 
produced by a distribution of matter-energy. This may 
be either a point-like particle, or an extended object. In 
both cases we assume that we are given the radial energy- 
momentum source function T^™"(r) for each mode [this 
is the source term in the s = +2 radial Teukolsky equa- 
tion (|^)]. For simplicity we assume here that the matter 
source is restricted to the range rmin ^ r < rmax (but 
this assumption may be relaxed - at least partially - as 
we discuss in section ^). 

Then i?^^'^(r) is given by 



2 '-'+2 \'^)^ 



R^_T{r) 



Tlr{r')H{ry)dr' , 



(69) 
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where 



H+{r,r') = a+{r')Rt2{r) 



H{r, r') = H+ (r, r')0{r - r') + (r, r')e{r' - r) , (70) 

and H^{r,r') are two smooth functions. We construct 
these functions from the two s = — 2 homogeneous ra- 



dial solutions R 



Xrauj{out) 



and ^ r: 



2' 



defined by their asymptotic behavior 

Rt2{r) o^r^e"^"^' (r - 

and 

i?:2(r) cx A2e-*'=" (r, - 



oo) 



-oo) ^ 



where k — oj ~ ma/ {2Mr^). (The s = +2 basis solutions 
are not required here. Also we do not require here a 
specific normalization for R^2 ^-^d RZ2-) We find 



H+{r,r')=a+{r')Rt2ir) 



and 



H- (r, r')^a- (r')i?l2 {r) + e-'^'' 



where u and r* are defined in Eqs. 
and 

a±(r') = {pW[R-_2.RU]y^ X 



(71) 

-*)^Br(,V, 
(72) 

|,El|), respectively, 



^(^')^^^2(^') + 5(r')i?f2(^') 



(73) 



Here p is a parameter given in Eq. (|2l|), and Ty[i?l2; -Ri2] 
is the Wronskian of the two basis functions (which is 
proportional to A), evaluated at r' . The functions 
A(r'), B(r') are specified in Appendix B. The four func- 
tions (r') are given by 



-a-(r') 



f,ir')RUr')+9.ir')^RUr') 
f,ir')RZ2ir')+g,ir')^RZ2{r') 



where fi and gi are functions of r' specified in Appendix 
A. 

The above construction yields the radial functions 
R^'^ir) for the solution that is valid (and regu- 
lar) through This domain includes the entire range 
T > Tmax, but not all points of r < rmax- The other 
solution ^~ that is valid through (i.e. everywhere 
in r < Tmin but not at all points of r > Tmin) may be 
constructed in a fully analogous manner. The only dif- 
ference is in the functions H^{r, r'), which now take the 
forms 



i=0 



(*- 



IX. SUMMARY OF MAIN RESULTS: 
ELECTROMAGNETIC PERTURBATIONS 

The electromagnetic case is treated in full analogy with 
the gravitational case. Here, again, the four-potential Aa 
is constructed in Chrzanowski's method, by applying a 
certain differential operator H^^ to a potential "^em- 



(74) 



Throughout this section we shall denote the electromag- 
netic potential '^em as ^I^ for brevity. This potential 
satisfies equations analogous to Eqs. (nn) and (O): 



and 



W-i [*] = 



(75) 



(76) 



where ifo is the s = +1 Weyl scalar, and W-i is the 
s = — 1 case of the differential operator (||). [See e.g. 
in which 5* is denoted 'V-e"- The last equation is the 
reduction of Eq. (15) therein to the Kerr case.] We use 
the decomposition 



and our goal is to construct the radial functions iZ^Y"- 
Again, we shall construct this function first in the case of 
pure electromagnetic waves, and then for perturbations 
with sources. Here we shall summarize the results; The 
main steps in the derivations are presented in Appendix 
C. 



A. Pure electromagnetic waves 

In this cases we assume that ipo is given. This field is 
decomposed into modes as 



J2 i?+r(r)5^T(6')e' 



For each mode Xmcu, R^'^ir) is a solution of the vacuum 
radial Teukolsky equation, and we assume this function is 
provided as a linear combination of two basis solutions. 
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The two sets of convenient basis solutions are: (i) the 
large-r set , 



and (ii) the ER set, 



Xmuidown) 



= Ae-''^" , ^ e''^" , (80) 



where k = to — ma/ (2Mr-^). 

Assume now that tp^ is given in terms of any two of 
the above four basis functions for i?^™'^, namely 

i?^T"(r) = A('^)i?^7"("^(r) + AWi?^7"^''^(r) , (81) 

and the coefficients A^"^ and A'^^^ are provided for each 
mode XmLO. (Here, again, "a" and denote two of 
the above four basis functions, e.g. ^'a''^="out" and 
" b''^ =" down" .) Then, the corresponding radial functions 
^Amw of ^ are given by 

i?^7-(r) = C('^)A('^)i?^'r('^^(r) + C(^)AWi?^™"(')(r) . 

(82) 

The four coefficients C^^"), C^™*), C'(''°"'"),C("p) take 
now the values 

C(™) = l/(4w^) , C(°"*) = 4tjVp (EM) , (83) 

C("P' = Q/p, C(''°"'") = 1/Q (EM) , (84) 
where in the electromagnetic case we have 

p = - 4a^w^ = + 4aw (m - aw) (EM) 

and 

Q^w{w + iq), (EM) (85) 

and, recall, 

w = AkMr+ , g = r+ - r_ = 2{M^ - c?fl'^ . 



B. Electromagnetic perturbations produced by 
sources 

Here we consider the case in which the perturbation is 
produced by charges and/or currents (e.g. a point charge 
or an extended charged object orbiting the BH). We as- 
sume that we are given the radial electromagnetic source 



function T^^(r) for each mode [this is the source term 
in the s = +1 analog of the radial Teukolsky equation 
(^]. As before we assume for simplicity that the source 
is restricted to the range Tmin <r < r^ax- 

The radial function R^^'^{r) then takes the form 



RtTir)= I T^T"{r')H{r,r')dr' , (86) 



where 



H{r,r') ^ H+{r,r')9{r - r') + H-{r,r')e{r' - r) , (87) 

and H^{r,r') are two smooth functions. We construct 
these functions from the two s = — 1 homogeneous ra- 
dial solutions = R+^ and = R-^^ 

defined by their asymptotic behavior 

Rtiir) oc re^'^''- (r -> oo) 

and 

i?;i(r) oc Ae-^'=" (r, ^ -oo) . 

(Here, again, we do not require a specific normalization 
for i?^i). We find 



H+{r,r') ^a+{r')Rtiir) 



(88) 



and 



H-{r,r') = a-(r')^li(0 +e"'^'""""'^*'Xl^»^(^')'' 



4=0 



(89) 



where u and r* are defined in Eqs. ( |42|j41| ), respectively, 
and 



a^{r') = -{pW[RZ„R+^]) ' 



A{r')—Rt,{r') + B{r')R^i{r') 



(90) 



Here W[R_i,R_^ = const is the Wronskian of the two 
basis functions R^i, and 

A{r') = 2iK , B{r') = A + 2iujr' - j l\ (EM) 

[with all quantities evaluated at r', e.g. K — am — (r'^ + 
a^)]. The two functions B~{r') are given by 



-a-(r') 



f,{r')R-_,{r')+g.{r')^RZ,{r') 



(for i = 0, 1), where the functions fi{r'), gi{r') are 
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/o(r') = (1 - iKr'/^) e*(™«-"'-*) , 

/i(r') = (VA)ife*('""-""'*) , 5i(r') = e*^""'-'^''*) 

(again, with u, r*, A all evaluated at r'). 

The above construction yields the radial functions 
R^'^(r) for the solution that is valid (and regu- 
lar) throughout S+. The other solution that is valid 
throughout E~ is constructed in a fully analogous man- 
ner. The only difference is in the functions H^{r,r'), 
which now take the forms 



H+{r,r') = a+(r')i?Ii(r) 



-i{mu—ujr-. 



(*- 



1=0 



H-{r,r') = a-{r')RZi{r) 



X. DISCUSSION 

Although in most of this paper we referred explicitly 
to gravitational perturbations, the same construction ap- 
plies to the electromagnetic case as well, as outlined in 
section [X. In particular, the domains of validity are the 
same in both cases: E"*" for (and for h'^^ or ^+ de- 
rived from the latter), and E" for (and for /i^^ or 

Also, although we have explicitly considered the ingo- 
ing radiation gauge throughout this paper, an analogous 
construction may be applied to the outgoing radiation 
gauge. In this latter gauge, too, there are two solutions, 
*Oi?G ^'^d *OflG (s-'^d correspondingly /iJ^q, ^qj^g 



which are valid in the domains E^^^ 



and Eq^p (but invalid in Eq^^, or Eq^q), respectively. 
The two domains Eq^^q are completely analogous to 
= Yifj^Q^ except that they are defined with respect 
to the ingoing rather than outgoing principal null con- 
gruence. In the rest of this discussion, too, we shall refer 
explicitly to the ingoing gauge, but the same remarks will 
be applicable to the outgoing gauge as well. 

Consider the case of a point particle. Our analy- 
sis shows there does not exist a single solution for the 
radiation-gauge hap or Aa that is regular in the entire 
off-worldline neighborhood of the particle. Instead, the 
solution '5+ (and correspondingly ft.^^,A+) has a line 
singularity along the outgoing null geodesic ^ emanat- 
ing from the particle towards the past and smaller r. 
Similarly, ^E*^ (and correspondingly h~^p^A~) has a line 
singularity along the null geodesic ^ emanating from the 
particle towards the future and larger r. The inevitabil- 
ity of such a line singularity in the radiation-gauge MP 



was previously demonstrated in Ref. [|j based on inde- 
pendent arguments. (The existence of ingoing radiation- 
gauge solutions other than /i^^,^^, which admit a line 
singularity in a different direction, not tangent to ^, has 
not been explored yet.) 

The unavoidable occurrence of a line singularity in the 
radiation-gauge fields ha^, A^ is obviously an inconve- 
nient property. Nevertheless it does not pose a too se- 
rious obstacle (at least in some important applications). 
We must recall that this singularity is after all a gauge 
artifact, which may in principle be removed by an appro- 
priate gauge transformation. Therefore, whenever the 
local values of hap or Aa are required for the calculation 
of any local gauge-invariant quantity, the solutions /i^^, 
yl+ and/or h~p, Aa may be used regardless of the fine 
singularity. 

An important application which requires the knowl- 
edge of hap or Aa is the radiation-reaction problem for 
a point mass or point charge. Generically the full anal- 
ysis of this phenomenon requires the calculation of the 
local self force acting on the particle. The electromag- 
netic self force is gauge-invariant. The situation in the 
gravitational problem is more delicate, because the grav- 
itational self force is a gauge-dependent entity. Neverthe- 
less, within the context of the adiabatic approximation, 
the orbit-integrated change (induced by the self force) in 
any of the orbit's constants of motion is gauge-invariant. 
One thus may use any gauge to calculate the self force, 
and hence the rate of change of the constants of motion. 
Consider the calculation of the self force according to the 
Mino-Sasaki-Tanaka jl^] formulation. Then the self force 
is the limit of the "tail- force" field at the particle's loca- 
tion. This limit may be taken from any desired direction. 
Two especially convenient directions are the ingoing and 
outgoing radial directions (so far the mode-sum method 
|]l3| has been fully developed for these radial directions 
only.) To this end, one may use the solution or A+ 
when calculating the self force from the radial direction 
r > ^particle, and the solution h~p or A'^ for calculating 

the self force from r < rparticia .In both cases the line 

singularity is not encountered. 



m 



In the case of a smooth extended source, E+ (or E^) 
becomes a four-dimensional set. In this case Vl'"'" does 
not develope an irregularity at E+; however, the equation 



***Recall, however, that in the gravitational case there is 
another difficulty associated with the radiation gauge: The 
leading-order asymptotic behavior of the MP, on approaching 
the particle's location from a generic direction, differs from 
that of the harmonic-gauge MP, jnaking this an "irregular 
gauge" in the terminology of Ref. |8| . This kind of irregularity 
(which is unrelated to the line singularity on E^) also occurs 
in e.g. the Regge- Wheeler gauge in the Schwarzschild case. 
This difficulty may in principle be overcome by transforming 
to an " intermediate gauge" , as outlined in Ref. |8| . 
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( p^ is violated there. This suggests that the quantity h'^^ 
(constructed from by applying the differential oper- 
ator n) will not be valid at S"*", even in its vacuum part 
(namely, it will fail to satisfy the vacuum Einstein equa- 
tion, and/or to reproduce the original Teukolsky field 
■00); But this still needs be verified. 

In the above construction we have assumed that the 
particle's worldline or the extended source is restricted 
to a range fmin ^ r < rmax- This assumption was made 
primarily for conceptual clarity, as it allows us to discuss 
the behavior of e.g. in the two vacuum regions, r > 
Tniax and r < r^in] but it can be relaxed at least partially, 
as we now discuss. 

Consider, first, the situation in which the source is re- 
stricted—to the range r < rmax with no minimal value 
Tniin. Then the construction of 4'"'" follows just as 
prescribed above, without any difficulties. The construc- 
tion of 4"" in this case may formally be carried out as 
above; However, the proof given in section [Vl| for the va- 
lidity of throughout E~ fails in this case: This proof 
(when applied to VP" rather than ^'^) starts from the 
trivial observation that (provided that the source is re- 
stricted to r > r„iin) Eq- (|l0| ) is satisfied hy ^~ through- 
out r < Tmiii. Then this feature is analytically extended 
to the entire domain T,~ . In the present case (i.e. no 
fmin) this proof is inapplicable even at its starting point. 
It therefore still needs be verified whether in this case the 
so constructed solution is valid in S~. 

In the analogous case, in which the source extends from 
infinity to some rmin, the situation is basically similar, 
though technically it is slightly more involved. Consider 
for example an unbounded orbit that arrives from infin- 
ity and scatters off the BH back to infinity. Here, the 
solution ^~ can in principle be constructed as above, 
but the solution is not guaranteed to hold (for the 
reason explained just above). In this case, however, due 
to the slow decay at large r of the potential term in the 
radial Teukolsky equation, the standard integral solution 
( p8| , p9| ) for -00 diverges. One then has to use another 
Green's function jl^ for the construction of ^Oj and this 
may modify the function H{r, r'). We shall not elaborate 
on this case here. 

Finally we note that there are a few types of spe- 
cial modes which require special treatment. First, for 
the stationary modes w = 0, the large-r basis solu- 



Xrauj{in,out) 
±2 



constructed in section IV must be re- 



tions R 

placed by some other ones, and the same holds for the 



^^^In the case of a point particle, this situation may be real- 
ized by a "fine-tuned" geodesic that asymptotes to an unsta- 
ble circular orbit in the far past, but falls into the BH in the 
future. We prefer not to consider here bound geodesies emerg- 
ing out of the white-hole horizon, to avoid the conceptual 
complications associated with the latter's causal properties. 



corresponding constants C'*^*"-' and C^°'^*\ Second, for 
" marginally-superradiant" modes fc = 0, the EH basis 
solutions /t_(_2 and the correspondmg constants 

(j{up,down) g^j,g |-|g modified. It appears likely, though, 
that in both cases the inhomogeneous solution described 



in e.g. section VIII remains valid, provided that one sub- 
stitutes the appropriate basis functions i?^2 ('") (i-C- those 
satisfying the correct boundary conditions at large r or 
at the EH). Other cases which require special attention 
are the so-called '7 = 0, 1 modes" (the '7 = mode" 
in the electromagnetic case) . These are the perturbation 
modes for which the Teukolsky variables "00 and 04 van- 
ish identically, while "02 is nonvanishing. The extension 
of this construction to include the '7 = 0, 1 modes", as 
well as all other a; = modes, is now underway. 
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APPENDIX A 



We rewrite Eq. ( pi) ) as 



(I?mc.)"[e-^(™-"^*)^B-(r')r'] 

1=0 

C+A+{r')(Dm.nR\{r)] 
-C-A-(r')(An.)"[i?l2(r)] 



(Al) 



(applied at r = r' and for n — 0, ...,3). It will 
be convenient to rewrite the polynomial TiiB~ {r')r'^ as 



Y.,B,{r'){r - r')\ By virtue of Eq. (|45D, the left-hand 
side of the last equation reads 



-i(mu—(jjr^ ) 



^i3,(r') j^(r-r')^ = e-'("™-"'--)n!B„(r') . 



1=0 



dr 



Evaluating Eq. (Al) at r = r' then implies (for n — 
0,...,3) 

e-'^("™""'--)n!B„(r') = C+A+(r')(D™^)"[i?l2] 

-C-A-{r'){D^^r[RZ2] (A2) 

(where D^^^^lR^^] is to be evaluated at r — r'). The 
operator D„iu is given by 



dr 
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Using the radial Teukolsky equation ( ]13| ) we can express 
any derivative ol an s = — 2 vacuum radial function i?^™'^ 
as a linear combination of i?^™'^ and {d/dr)R^2'^ (and 
consequently we can express any power of Dmuj as a lin- 
ear combination of i?^^'^ and -Dmc^[-R^™"])- As a con- 
sequence, when applied to any homogeneous solutions 
(and in particular i?^2)i have the following op- 
erator identities: iisl 



[D^u^f = {2/A){iK + r - M)D„,^ + 



d 



A + 6iujr 



(2/A)(iK + r - M) — 
dr 

+A''^[2iK{iK + r~M) + A(A + 6iLur)] 



and 



{D,nujf = A-^[4:iK{iK + r - M) + (A + 2 - 2iujr)A]D,„ 
+A'~^[2iK{\ + Giujr) + 6iujA] 

= UiKiiK + M) + (\ + 2- 2iujr)A] 

dr 

+A-^{iK[4:iK{iK + r- M) 
+ (A + 2-2iwr)A] 
+A[2iA'(A + Qiujr) + QiujA] } . 

Thus, we may write Eq. (|A^) as 



Bn{r') = C+A+{r') 
-C-A-{r') 



(for n = 0, 3), where the functions /„, g„ are given by 



, ga{r) ^ 0, 



/i(r) = (i/A)Ke*(""'-"'^*) , 51 (r) = e*'""""'^*) , 



/2(r) = (2A2)-i[2iX(iX+r-M)+A(A+6iwr)]e*("™-'^''-) , 



52(?') = A ^(iif + r-M)e 



/3(r) = {^A^y^iiKyiiKiiK + r - M) 
+ (A + 2-2iwr)A] 

+A[2?;if(A + Uujt) -f 6it^A] }e^("™-'^'-*) , 

g-iir) = {^A^)-^yiiK[iK + r - M) 

+ {X + 2-2iu;r)A]e'^"'''-'^'''K 

Once the coefficients Bi{r') are determined, the origi- 
nal coefficients B~{r') may be constructed through 



^E-(r'V = ^S,(r')(r-r'r, 



1=0 

which yields 



i=0 



= Bo - Bir' + B^r'^ - B3r'\ 
=Bi- 2B2r' + 3S3r'2 , 
B,- = ^2 - SSar' , B:^ ^ 4 . 



f,,{r')RUr')+gn{r')^RUr') 



This allows us to express the functions B^ (r') as 
B-ir') = C+A+ir') 

-C-A-{r') Ur')R-_2{r')+gn{r')^RZ2{r') 

with the functions fnij"'), gnif') given by 

/o = /o - fir' + hr'^ - /3r'3 , /i = A - 2/2r' + S/gr'^ , 
/2 = /2 ^ S/ar' , fs — fs : 

and similarly 

50 = 50 - .gi''' + 52?"'^ - 53?^'^ , 51 = 51 - 2g2f ' + Sgar'^ , 
52 = 52 - 3.g3r' , 53 = 53 

(with all functions gn, /„ evaluated at r' rather than r). 

APPENDIX B 

Our goal is to construct the basis functions i?^2 from 
the corresponding functions R^2- To simplify the nota- 
tion, throughout this appendix we shall view i?^2 and all 
other "radial" variables as functions or r, not r' . When 



implementing the result (B6) back in section VI, one 
should simply substitute r r' . 

The analysis in section [V, Eqs. ( p7| , p2[ ), implies 



which [since H is the inverse of the operator (I?mtj)^] is 
equivalent to 

R+2 ~ {DjnLj)'^[R^2\ ■ 

With the aid of the radial Teukolsky equation, the op- 
erator {D„u^)* acting on any vacuum solution R'^'^ 
may be expressed in terms of R'^'^ and its first-order 
derivative. Chandrasekhar derived the formula [see Eq. 
(49CH) , where throughout this Appendix " CH" refers to 
equations in chapter 9 of Ref. ||]] 

(A„c.)"[i?^'r] = (Ao/A3)An^[i?^r] + (So/A3)i?^r , 

(Bl) 



18 



where 



APPENDIX C 



Aq = -BiKlK"^ + (r - M f] 

+ [AiK{\ + 2) - 8itjr(r - m)]A + i 



(B2) 



Bo = [(A + 2 - 2iwr)(A + Qiujr) + UicoiiK -r + M)]A 
+AiK{iK - r + M) (A + Qiuor) . (B3) 

This yields 

i?±2 = (C7±/A3) (AoI?™„[i?±2] + B^R%) . (B4) 

We also need to express the determinant Ty[i?^2J -^+2] 
in terms of R^2- We find it useful to express this deter- 
minant as 

W[R+2i ^+2] = -R+2-R+2,r - -R+2-R+2,r 

= -R+2^-l[-R+2,r] - ^+2,r^-l[-R+2] ' 

where 2?Li = dr - (iK/A) - 2{r - M)/A, which allows 
us to make use of Eq. (50CH). Writing 

vl^iAoD^^ + Bo) = AiD^^ + Bi 

[with Ai,Bi specified in Eq. (51CH)], we obtain 

W[R-2,Rt2] = C+C- A-^{BaAi- AoBi)W[RZ2, Rt2] ■ 
A straightforward calculation yields 

BoAi - AoBi =. pA^ , 

leading to 



W[R^2,K2] = C+C-pA-^W[RZ2,R±2 



(B5) 



[Note the consistency of this result with the general ex- 
pression for the Wronskian of the Teukolsky equation: 
For any s, and any pair of independent solutions i?^, 



T^[i?" i?^] ^ const -A 



-s-l 



Hence W^[-R^2:^+2] = const ■ A~^ and W[RZ2,R-2\ = 
const ■ A, in agreement with Eq. (B5).] Combining this 
result with Eqs. (B4) and (p3), we obtain 



A± = {C^pW[RZ2,Rt2]) {AoDm^[RT2] + BoRZ 
This yields 



BR^. 



with 



A = Ao , B = Bo + {iK/A)Aa . 



(B6) 
(B7) 



For a vacuum mode \muj of electromagnetic pertur- 
bations, the radial function i?^™" must satisfy the two 
equations 



pAmc^j^Amo^^^^j =0 



and 



(CI) 
(C2) 
(C3) 



R+r{r) = -{D,n^?[R^^^{r)]. 
The general solution for these two equations is 

RtT ^ -p-'A{Dl^fA [i?ir (r)] , 
where 

p = A^ - Aa^uj"^ = A^ + Aauj (m - auj) (EM) . 



Considering the four asymptotic basis solutions 
/c_j_2 ^ ' ' specmed m section |LXj, the corre- 



sponding parameters C("\ C(°"*), C^down) ^(jiup) are eas- 
ily calculated just as in the gravitational case. One finds 

= 1/(4^2) , C(°"*) = iuj^/p (EM) 

C7("^'' - Q/p , C(''°"'"' = 1/Q (EM) . 

The general solution to the homogeneous part of Eq. 
, namely {Dmui)'^[R^^] = 0, is easily constructed: 



^^T"(r) = e-*(™"-'^''*) Bi 



(C4) 



Consider next the case of inhomogeneous electromag- 
netic perturbations. The general solution for the radial 
function of tpo may be expressed as 



(C5) 



The Green's function is 

G{r,r') = A+{r')Rlj^{r)9{r~r')+A-{r')R-j^{r)0ir'-r) , 
where 

A^=RIJ{AW[R^„R+,]) (C6) 
(evaluated at r'). Then R^^{r) is given by 

R^^'^{r)^ J T^™'^(/)B(r,r')dr', (C7) 
^"111111 

where H{r,r') satisfies 



19 



{Dr^u f [H{r,r')] = -G{r,r') 



(C8) 



(and the appropriate boundary conditions at r > rmax)- 
We find H{r, r') to be of the form 

H{r,r') = H+{r,r')9{r - r') + H-{r,r')e{r' - r) , (C9) 
with 

H+{r,r') = C+A+{r')R\ir) (CIO) 

and 

1 

H- (r, r') = C-A- (r')i?Ii (r) + e"*'™"-"'^-) ^ Br (r'y 

i=0 

(Cll) 

The two functions B~ {r') are determined by regularity 
conditions at r = r\ which yield 



Brir') = C+A+ir') 
-C-A-{r') 



f,{r')RU{r')+g,{r')^RU{r') 



f,{r')R-_,{r') + CH{r')^R-_,{r') 



In full analogy with the gravitational case (see Appendix 
A) we find 

fo = fo- fir' , /i = /i , 50 = 50 - 9ir' , 5i = 5i , 
and 

/o(r) = e^(™™*) , 5o(r)=0, 
/i(r) = (i/A)Ke*(™"-"'^*) , 51 (r) = e'(™"-'^'-*) , 



yielding the functions fi , gt specified in section [X 
Next we express R^^ in terms of R^i, using 



R^i — — C {D„iuj) 



(C12) 



The vacuum radial Teukolsky equation for the s = — 1 
radial function may be expressed as ^ 

Dl^D„,^ = (A + 2iujr)/A . 

Using this equation to reduce the order of differentiation, 
and recalling Dmuj = D}^^ + 2iK/ IS., we obtain 

{D.mu:f[R^l] = (^o/A)A„c.[i?^i] + (So/A)i?±i , 

(C13) 



i?^i = -(C±/A)(A-i?±,+Bi?±J, 



where 



and 



2iK 



(C14) 
(C15) 



B = Bo + {iK/A)Ao = A + 2iujr ~ 2K^/A . 
We now calculate the determinant W^[i?^]^, i?^]^], using 

W[R^i, R'^i] = R^iDm^[R'^i] — R'^iDmuj[R+i] , 
along with Eqs. ( C12 , cT^ ). The calculation yields 
= C+C-pA-^W[RZi,Rti] . 



Substituting this and Eq. (C14) into Eq. ( |C6|) we obtain 
a±(r) = C±^±(r) 

= -{pW[RZ„Rt,])-\A-^Rf, + BRf,) . 

Finally, substituting this in the above equations for 
B^{r') and H^{r,r') (with the substitution r — » r'), we 
obtain the expressions for these quantities as specified in 
section 
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